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ABSTRACT 

We have performed integral field spectroscopy of the planetary nebulae Hen 3-1333 
(PNG332.9—09.9) and Hen 2-113 (PNG321.0+03.9), which are unusual in exhibiting dual¬ 
dust chemistry and multipolar lobes but also ionized by late-type [WC 10] central stars. The 
spatially resolved velocity distributions of the Hct emission line were used to determine their 
primary orientations. The integrated Ha emission profiles indicate that Hen 3-1333 and Hen 2- 
113 expand with velocities of ~ 32 and 23 km s _1 , respectively. The Hubble Space Telescope 
observations suggest that these planetary nebulae have two pairs of tenuous lobes extending 
upwardly from their bright compact cores. From three-dimensional geometric models, the 
primary lobes of Hen 3-1333 and Hen 2-113 were found to have inclination angles of about 
—30° and 40° relative to the line of sight, and position angles of —15° and 65° measured east 
of north in the equatorial coordinate system, respectively. 

Key words: stars: Wolf-Rayet - ISM: kinematics and dynamics - planetary nebulae: general. 


1 INTRODUCTION 

The planetary nebulae (PNe) Hen 3-1333 (= PNG332.9—09.9) and 
Hen 2-113 ( = PNG321.0+03.9) show unique properties of the 
early post-asymptotic giant branch (AGB) phase and dust for¬ 
mation at the last AGB phase. ICohen et afl i 19991. l2002h found 
a dual-dust chemistry in these PNe using infrared observations. 
ICohen et all ( 1999 1 fitted blackbody curves to the continuum of 
Hen 3-1333, and found some evidences for dual-chemistry, namely 
carbon-rich and oxygen-rich dust grains. Similarly, ICohen etaD 
( 2002 1 identified oxygen- and carbon-rich materials in Hen 2-113. 

The central stars of Hen 3-1333 and Hen 2-113 were among 
the coolest Wolf -Rayet stars studied in JHKL bands and classi¬ 
fied as [WC 11] (I Webster & Glassl IT974h . They have been classi¬ 
fied as lat e-type [WC10] based on the classifi c ation scheme pro¬ 
posed by [Crowther et 5] dl998h . [Gleizes et afl dl989h derived the 
H I Zanstra temperature of T z (Hl) = 17.5 kK for Hen 3-1333 
and 37kK for Hen 2-113. Furthermore, photoionization modelling 
bv iDe Marco & CrowtheJ dl998l) yielded T e ff — 25 kK for Hen 3- 
1333 and T eff = 29 kK for Hen 2-113. 

FigJU shows the Hubble Space Telescope (HST) images of 
Hen 3-1333 and Hen 2-113 taken with the High Resolution Chan¬ 
nel of the Advanced Camera for Surveys (ACS/HRC), and through 


* Based on observations made with the Australian National Univer¬ 
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the F606W and T814W filters, respect i vely. T he HST image of 
Hen 3-1333 studied bv IChesneau et alJ d2006h hints at complex 
mul tipolar lobes surrounding its bright compact core. Previously, 
IDe Marco et alJ ( 20021) identified a c ompact dusty disk in Hen 3- 
1333. Moreover, Lagadec et alJ d2006h described the HST image of 
Hen 3-113 as two ring-like structures produced by the projection of 
a hourglass-shaped geometric model. 

In this paper, we present spatially resolved kinematic observa¬ 
tions of Hen 3-1333 and Hen 2-113 made with an integral field unit 
(IFU) spectrograph. We first describe the observational method and 
the data obtained, and then proceed to determine the spatial orien¬ 
tations constrained by kinematic models. 


2 OBSERVATIONS 

The IFU observati ons were performed using the Wide Field Spec¬ 
trograph (WiFeS: Ibopita et aD 120071. [201(1) mounted on the 2.3- 
m Australian National University (ANU) telescope on 2010 April 
20. TableQ]presents an observation journal, including the exposure 
time used for each PN in our WiFeS observations (column 3), and 
information on the HST observations (columns 5-8). 

WiFeS is an image-slicing IFU developed and built for the 
ANU, feeding a double-beam spectrograph. It has a field-of-view 
(FOV) of 25 arcsec X 38 arcsec and a spatial resolution of 1 arc- 
sec. We used the spectral resolution of R ~ 7000. The classical 
data accumulation mode was used, so a suitable sky window has 
been selected from the science data for the sky subtraction pur- 
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Table 1. Journal of observations. 


Object 

Other name 

WiFeS 



HST 




Exp. (s) Obs. Date 

Filter 

Exp. (s) 

Obs. Date 

Programme ID 

Hen 3-1333 

PNG 332.9-09.9 

1200 2010 Apr 20 

F606W 

56 

2002 Sep 17 

9463 

Hen 2-113 

PNG 321.0+03.9 

60,1200 2010 Apr 20 

F814W 

56 

2003 Mar 09 

9463 


Hen 3-1333 


S lobe 


Hen 2-113 


NE lobe 


N lobe 



SW lobe 


NW lobe 



1 ” 


Figure 1. HST images of Hen 3-1333 (top panel) taken with the F606W 
filter, and Hen2-113 (bottom panel) taken with the F8\4W filter and 
ACS/HRC detector (Programme ID 9463, PI. R. Sahai). The image scale 
is shown by a solid line in each image. North is up and east is towards the 
left-hand side. 


pose. We also acquired series of bias, dome flat-field frames, twi¬ 
light sky flats, arc lamp exposures, wire frames and standard stars 
for bi as-su btraction, flat-fielding and calibrations (fully described 
by|D ane hkar2014). 

The WiFeS data were reduced using the IRAF pipeline WIFES 
(version 2.14). We ut ilized the same d ata red uctio n procedure de¬ 
scribe d in detail by iDanehkar et alJ d2013h and iDanehkar et alJ 
(120141) . which includes flat-fielding, bias-subtraction, wavelength 
calibration, spatial calibration and flux calibration. 
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Figure 2. From left to right, spatial distribution maps of flux intensity and 
LSR velocity of Ha A6563 for (a) Hen 3-1333 and (b) Hen 2-113. Flux 
unit is in 10“ 15 erg s -1 cm -2 spaxel -1 and velocities in kms -1 . The 
white/black contour lines show the distribution of the narrow-band emis¬ 
sion^ Ha in arbitrary unit obtained from the SuperCOSMOS Ha Survey 
^Parker et~afll20051) . North is up and east is towards the left-hand side. 


3 OBSERVATIONAL RESULTS 

Fig-El presents spatially resolved maps of flux intensity and radial 
velocity derived from the Ha A6563 emission line for Hen 3-1333 
and Hen 2-113. To extract them, we utilized a Gaussian curve fit¬ 
ting code, which allowed us to fit a single Gaussian profile to the 
Ha emission line for all spaxels across the IFU field. As seen in 
Fig. [2] the radial velocity maps indicate the presence of outer faint 
lobes with a diameter of about 10 arcsec extending to the compact 
structure of ~ 3.5 arcsec diameter (see Fig.Q]). Contour lines in the 
figure depict the 2P distribution of the Ha emissi on obtained from 
the SuperCOSMOS Ha Sky Survey ( Parker et al . 2005 ), which 
may aid us in distinguishing the nebular borders. However, the 
broader point spread function (PSF) of the stellar Ha emission cre¬ 
ates diffraction spikes around these PNe, so the contours do not 
actually show the nebular borders. The IFU flux map presents the 
resolved nebular Ha emission. 

Table [2]hsts our velocity results for different emission lines of 
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Figure 3. HST images at two different contrast levels (left-hand and second panels), SHAPE mesh model (third panel), and rendered image (righ-handt panel) 
of (a) Hen 3-1333 and (b) Hen 2-113. The mesh models are shown at the best-fitting inclination. The synthetic rendered images have been obtained using a 
Gaussian blur. 


Table 2. LSR systemic velocities and HWHM expansion velocities. 


Name 

Z’sys(Ha) 


kHWHM(kms 

- 1 ) 


(km s _1 ) 

Ha 

[N II] 

[S II] 

Mean 

Hen 3-1333 

-62.2 

31.6 

37.2 

30.6 

33.1 ±3.3 

Hen 2-113 

-56.7 

22.5 

22.3 

20.1 

21.6 ± 1.2 


the integrated spectrum of each PN. Column 2 presents the local 
standard of rest (LSR) systemic velocity ( v sys ) derived from the 
Ha emission line. Our expansion velocities (u exp ) derived from 
the half width at half-maximum (HWHM) for Ha A6563, [N II] 
AA6548,6584 and [S II] AA6716,6731, and average HWHM values 
are presented in columns 3-6, respectively. Our observations for 
the Ha emission line give Vhwhm = 31.6 and 22.5 km s -1 for 
Hen 3-1333 and Hen 3-113, in good agree ment with u exp — 30 
and 19 km s" 1 derived by De Marco et alJ dl997h . respectively. For 
Hen 3-113. lGesicki et al1 ( 20061) also derived an expansion velocity 
of 18 km s- 1 and a turbulence velocity of 15 km s -1 . 

From the HST image of Hen 3-1333 displayed in the top panel 
of Fig.[l] several elements can be identified: two pairs of bipolar 
lobes in the north and south, and some filamentary structures. The 
HST image of Hen 2-113 (shown in the bottom panel of Fig.Q]), 
shows two pairs of bipolar lobes in the east and west, and a dark 
lane between the western border of two opposing lobes. But, the 
HST images are also obscured by the PSF diffraction spikes, which 
need to be considered carefully when studying the morphology. 


4 THE MODELS 

The 3D kinematic mode lling program SHAPE dSteffen & Lopezl 
l200d:ISteffen et al.ll201 ih was used to identify the spatial distribu¬ 
tion of the outer tenuous lobes around the compact cores of Hen 3- 
1333 and Hen 2-113. The modelling procedure consists of defining 
a geometry, assigning emissivity distribution and defining a veloc¬ 
ity law as a function of position. The program produces several out¬ 
puts that can be directly compared with the observations, namely 
position-velocity diagram and appearance of the object on the sky. 

For this study, we adopted a bipolar geometric model (marked 
as a — a' in Fig.[3j left-hand panel). To replicate the HST images, 
another bipolar component, which is identical to the other compo¬ 
nent, is introduced with different orientations (marked as b — b' 
in Fig.0>. We adjusted the inclination (z) and position angle ‘PA’ 
in an iterative process until the qualitatively fitting solution is pro¬ 
duced. The velocity IFU maps derived from the Ha emission line, 
combined with the HST images, were used to constrain the primary 
bipolar lobes (a — a'). The HST observations were also used to re¬ 
construct the secondary bipolar lobes (b — b'). Fig. [3] (second panel) 
also shows the HST images with higher contrast, which reveal de¬ 
tails not accessible at the previous contrast level. 

Fig-0 (third panel) shows a 3D representation of the model 
at the best-fitting inclination of each object. From the best-fitting 
kinematic models, the primary lobes (a — a') of Hen 3-1333 and 
Hen 2-113 were found to have PA = —15° ± 5° and 65° d= 5°, 
and inclination angles of z = —30° d= 15° and 40° =b 15°, respec¬ 
tively. The secondary lobes (b — b') of Hen 3-1333 and Hen 2-113, 
which are only noticeable in the HST images, were found to have 
PA = 60° =b 5° and 120° zb 5°. Fig.0 (right-hand panel) shows the 
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Table 3. The key parameters of the best-fitting kinematic models. 




Hen 3-1333 



Hen 2-113 


Feature 

Size (arcsec) 

PA(°) 

GPA(°) 

i(°) 

Size (arcsec) 

PA(°) 

GPA(°) 

i(°) 

a — a' 

6.2 x 3.7 

-15 ±5 

39.6 ±5 

-30 ± 15 

3.7 x 2.4 

65 ±5 

93.3 ±5 

40 ± 15 

b-b' 

6.8 x 2.9 

60 ±5 

114.6 ±5 

_ 

4.8 x 2.6 

120 ±5 

148.3 ± 5 

_ 


synthetic rendered images. As seen, the rendered images perfectly 
resembles the HST observations shown in Fig. [3] (left-hand panel). 

Table [3] lists the key parameters of the best-fitting morpho- 
kinematic models: the sizes, the PA, the Galactic position angle 
(GPA), and the inclination ( i ) of the lobes, respectively. The PA 
is the position angle of the bipolar lobes projected on to the plane 
of the sky, and measured from the north towards the east in the 
equatorial coordinate system (ECS). The GPA is the position angle 
projected on to the sky plane, measured from the North Galactic 
Pole (NGP) towards the Galactic east. The inclination is measured 
between the line of sight and the nebular symmetry axis (i = 0° 
being pole-on). 


5 CONCLUSION 

The spatially resolved kinematic observations of Hen 3-1333 and 
Hen 2-113 presented in this paper have allowed us to identify their 
primary orientations. Our kinematic maps also indicate that these 
PNe have large extended faint lobes upwardly from their com¬ 
pact structures. For the primary lobes of Hen 3-1333 and Hen 2- 
113, we derived the inclination angles of i — —30° and 40°, re¬ 
spectively. Accordingly, the main orientations of Hen 3-1333 and 
Hen 2-113 were found to be PA = —15° and 65° in the ECS , 
in exc ellen t agreement w i th the HST studies by Che sneau et alJ 
(2006j) and lLagadec et alJ (l2006h . respectively. Interestingly, both 
Hen 3H333 and Hen 2-113 have the same stellar characteristics 
dDe Marco & Crowthej[l998 ). but they show different expansion 
velocities (see Table [2]). Therefore, it seems that their nebular kine¬ 
matic features are somehow unrelated to their stellar characteris¬ 
tics. 

Both Hen 3-1333 and Hen 2-113 demonstrate a dual-dust 
chemistry consisting of c arbon-rich and oxyg en -rich grains 
( Cohe n et aD 1 19991 . I 2 OO 2 I : iDe Marco et all l20oi) . ICohen et alJ 
(2002) also found the sam e properties i n oth er PNe with late- 
type [WC] stars. Moreover, iGornv et aD (2010) found more PNe 
with dual-dust chemistry in the Galactic bulge, and speculated 
that the simultaneous presence of O-rich and C-rich dust grains is 
more likely related to the stellar evolution in a close binary sys¬ 
tem. Recently, iGuzman-Ramirez et al.1 ( 2014) identified the pres¬ 
ence of dense central tori in PNe with dual-dust chemistry, sug¬ 
gesting the possible formation through a common-envelope phase. 
Meanwhile, some recent kinematic studies of bipolar PNe around 
close-binary central stars indicate that their binary orb ital incli na¬ 
tions are very close to their nebular inclinations (see e.g. l Jones et alJ 
l2012l : lTvndall et al.ll2012l : lHuckvale et aPl2013l) . Therefore, deeper 
observations of the central stars of Hen 3-1333 and Hen 2-113 will 
lead to a better understanding of their bipolar morphology and dual¬ 
dust chemistry. The presence of a stellar companion should be in¬ 
spected. However, it is extremely difficult to detect a substellar 
companion, which could also affect the nebular properties. 
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